Quantum dynamics in a camel-back potential of a dc SQUID 
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We investigate the quantum dynamics of a quadratic-quartic anharmonic oscillator formed by a 
potential well between two potential barriers. We realize this novel potential shape with a supercon- 
ducting circuit comprised of a loop interrupted by two Josephson junctions, with near-zero current 
bias and flux bias near half a flux quantum. We investigate escape out of the central well, which can 
occur via tunneling through either of the two barriers, and flnd good agreement with a generalized 
double-path macroscopic quantum tunneling theory. We also demonstrate that this system exhibits 
an "optimal line" in current and flux bias space along which the oscillator, which can be operated 
as a phase qubit, is insensitive to decoherence due to low- frequency current fluctuations. 



Superconducting devices, based on the nonlinearity 
of the Josephson junction (JJ), exhibit a wide variety 
quantum phenomena. During the last decade, inspired 
by Macroscopic Quantum Tunnelling (MQT) studies [Ij, 
quantum dynamics of the current biased J J, dc SQUID 
and the rf SQUID phase qubit have been extensively 
studied |2] [3l [4l [S] [6] . In each of these devices, the dy- 
namics can be described as those of a quantum particle 
in a quadratic-cubic potential. The flux qubit system [7], 
realized by three or four JJs in a loop, is described by a 
double well potential. Here we propose to study a new 
potential shape called hereafter a "camel-back" double 
barrier potential, shown in Fig. [T]^. This potential is ob- 
tained using the dc SQUID circuit shown in Fig. [l^ in 
a new way. The characteristics of the camel-back po- 
tential, including depth and relative barrier height, are 
controlled by the SQUID current bias and flux bias 
$ext- There is a special line we call the "optimal line" 
in these two bias parameters at which the barrier heights 
are equal and anharmonicity is quartic. Because of the 
symmetry of the potential, the system can escape from 
the central well via tunneling through either of the two 
barriers to an adjacent deeper well. We investigate this 
double path escape and analyse it using MQT theory. 
When the depth of the central well is adjusted so that 
the escape rate is negligible but anharmonicity is signif- 
icant, the two lowest energy levels |0) and |1) constitute 
a phase qubit. We show that this qubit is insensitive to 
dephasing due to current bias fluctuations on the optimal 
line. 

A dc SQUID circuit has two degrees of freedom corre- 
sponding to the phase differences <pi and (j)2 across its two 
JJs. The dynamics are analogous to those of a particle 
of mass m = 2C($o/27r)2 in the 2-D potential [HIS] 
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FIG. 1: Schematic of experimental setup and camel-back po- 
tential, (a) Circuit layout, (b) Full 2-D potential for b — 3.05, 
?7 - 0.72, a = 0, $ext - -0.508^0, /b - 0, showing the fam- 
ilies of minima associated with the [0$o] and [— l$o] fluxoid 
states. The black line follows the minimum energy path, (c) 
Potential along the minimum energy path, parameterized by 
the path length. 



for a given sample are the Josephson energy Uq = 
{Id + /c2)*^o/27r, the junction to loop inductance ra- 
tio h = <I>o/27rL/c, the critical current asymmetry a = 
(/c2 — ^ci)/2/c7 and the loop inductance asymmetry rj = 
(L2 — Li)/L. Here Ic = Id + ^c2, hi and /c2 are the 
critical currents of the two junctions, Li and L2 are the 
geometric inductances of the two arms of the SQUID 
loop, L = Li +L2, C is the capacitance of each junction, 
and $0 = h/2e is the quantum of flux. The external 
control parameters and <I>ext enter into the potential 
through yi) = 7r<I>ext/^o and s = h/Ic- For our sample, 
Ic = 11.22 /iA, C = 250.3 £F, 6 = 3.05, r] = 0.72, and 
a = 0.0072. 

The first term in U{x^y)^ due to the junctions, de- 
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scribes a 2-D periodic array of minima and maxima. This 
array can be tilted in the x-direction with an apphed 
current bias. Magnetic field energy associated with cir- 
culating current gives rise to the parabolic term in the 
^/-direction, the minimum of which is shifted by the ex- 
ternal fiux. 

Stable, stationary states of the system correspond to 
minima of U{x^y). There can exist one, two, or more 
minima families corresponding to distinct fiuxoid states 
[n<l>o]. For each, when s exceeds a fiux dependent critical 
value 5c[n<l>o](^6), the related local minima disappear. 
For small values of 6, the parabolic term in U (x, y) is 
shallow, and there can be many fiuxoid states. For h ^ 
I/tt, as in our case, the parabolic term is steep and there 
is only one stable fiuxoid state except in a small region 
around <l>ext/*^o = 0-5 (mod 1), ~ where there are 
two states with opposite circulating current. Hereafter 
we will be focusing on this region. 

In general, dynamics is described by 2-D motion in the 
potential. In our case, the particle moves through a val- 
ley in which the curvature is much larger in one direction 
(-100 GHz) than the other (10-20 GHz). To a good ap- 
proximation therefore, motion is one dimensional along 
the path of minimum curvature which connects minima 
and saddle points. We parametrize this path with the 
phase length z (see black line in Fig.[T]3). U{z) in Fig.jlJ^ 
depicts the "camel-back" potential shape we are investi- 
gating. In a typical experiment, the system is initialized 
in the central well, which corresponds to the [O^o] fiux- 
oid state. The deeper wells on either side of the central 
well both correspond to the [-l<^>o] fiuxoid state. Starting 
from the central well, the system can escape via tunnel- 
ing through the barriers in either of the two physically 
distinct directions to the [-I^q] fiuxoid state. 

In the perfectly symmetric case, the potential near the 
central minimum will be harmonic with a quartic pertur- 
bation. More generally, the Hamiltonian for small oscil- 
lations in U {z) is 

H = ^hujp{P^ + Z^) - ahujpZ^ - dhujpZ^. (2) 

Here cOp is the zero amplitude oscillation frequency in the 
direction of minimum curvature, and Z = z^mujpjli and 
P = p/ ^yh^Jpm are the reduced position and correspond- 
ing momentum operators. Treating the anharmonic 
terms as perturbations, to second order the transition en- 
ergy between levels n-1 and Tl is flLUfi — i^fi — huj(^l — 72 A) , 
where the anharmonicity is A = + 3(5 [10 . We 

have calculated the escape probability for the camel- 
back potential with a double escape path in the quantum 
limit using the instanton formalism [11 . For a dura- 
tion At, it reads Pescjh ^ ^ext) = 1 - e-(^^+^^)^^ where 
^R,L = Ar^l^^/Nr^l exp[-BR^LNR^L]' Here R and L 
refer to the right and left barriers. Nr^l = ^Ur^l/^^ are 
the normalized barrier heights. The general expression 
of the coefficients Ar^l, and Br^l depends on the poten- 



tial shape. In the symmetric case where cr(/5, ^ext) = 0, 
the potential is quadratic-quartic, Ar^l = 22 7r~2 and 
Br^l = 16/3. Far from this symmetric line the potential 
is quadratic-cubic, the escape rate through one barrier 
is dominant {e.g. = 0), and we retrieve the standard 
MQT situation {S = 0): Ar = 6^71'^ and Br = 36/5 [J. 
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FIG. 2: Escape from the ground state, (a) Critical lines 
of three neighboring fiuxoid states, measured (symbols), and 
standard MQT theory fit (lines), denoting /5o%, the ampli- 
tude of a 60 /is h pulse that yields Peso = 50% to the volt- 
age state of the SQUID, (b) Critical lines representing 50% 
escape out of fiuxoid states [0$o] and [-l$o] in the region 
around $ext = — $o/2 due to a 100 ns flux pulse, measured 
(symbols), and the generalized MQT theory flt (lines). The 
[0$o] and [-l^o] fluxoid states are both stable in the central 
region enclosed by the critical lines. Escape is either to the 
SQUID voltage state {h > 0.8 /xA) or to the adjacent fluxoid 
state {h < 0.8 /xA). The points W, M, and R indicate the 
Working point, quantum Measurement point and Readout 
point for a typical camel-back potential phase qubit exper- 
iment, (c) Width of ground- state escape A$: measurements 
(points-h lines) , and generalized MQT theory with (solid line) 
and without (dashed line) 9 nA RMS low-frequency current 
noise. The location of the dip near the maximum A$ corre- 
sponds to the point where symmetry leads to a reduction in 
sensitivity to noise. 

A schematic of our experimental setup is shown in 
Fig. [l^. Our sample was fabricated at PTB using a 
Nb/AlO^c/Nb trilayer process with Si02 dielectric and 
a critical current density of 300 A/cm^ [12]. The 5 /im^ 
junctions are embedded in a square loop with inner size 
10x10 /im^. An off-chip coil provides a dc fiux bias. 
Current bias and voltage leads, heavily filtered at var- 
ious stages of the cryostat[13^, connect at the right of the 
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SQUID. Fast flux pulses are inductively coupled via the 
on-chip loop to the left of the SQUID. Microwave (MW) 
excitation is applied via an on-chip loop which couples 
inductively to the current bias leads. The MW excita- 
tion must be in current, rather than flux, because for the 
symmetric camel potential, small amplitude oscillations 
occur for the most part in the x direction, and therefore 
must be excited via the —sx term in U{x,y). The fast 
flux and MW excitation lines are 50 Q coaxial with -20 dB 
attenuators at 1 K and base temperature. The SQUID 
chip is enclosed in a small copper box thermally anchored 
to the mixing chamber of a dilution refrigerator with a 
base temperature of 30 mK. The cryostat is surrounded 
by a superconducting Pb shield, inside a /i-metal shield, 
inside a soft iron shield. 

Fig. shows the switching current c:^ 

IcSc[n^o]{y})) as a function of flux for the [-l<^>o], [O^o], 
and [l^o] fluxoid states. The interior of each curve is 
the region where the corresponding flux state is stable. 
The measurements shown in Fig. [2^ were obtained with 
a standard technique in which pulses of varying am- 
plitude are applied and a dc voltage detected across the 
SQUID when it switches to its voltage state. With this 
scheme, however, there is no direct indication of multiply 
stable flux states. In Fig. we use a novel technique to 
measure the overlapping critical lines of [O^o] and [-l<^>o] 
flux states close to $ext/^o = —0.5. These two interior 
critical lines represent transitions between the two flux 
states, rather than transitions to the voltage state, which 
is why the standard technique does not detect them. 

Our novel escape measurement method proceeds as fol- 
lows. First, if necessary, the system is initialized in the 
desired flux state with an adiabatic pulse on the fast flux 
line. I}) is brought to its working point value. A flux pulse 
is applied via the fast line for a flxed nanosecond-scale 
duration, bringing the total externally applied flux to a 
"measurement point" close to the critical line. This has 
the effect of reducing the heights AUr^l of the two po- 
tential barriers. Pesc via tunneling from the central well 
through the barriers to the neighboring deeper wells is 
thereby greatly increased. The system is brought back 
to a flux at which both fluxoid states are stable. The 
fluxoid state is then read out via a slow (~ 10 /is) 
pulse. This pulse brings the system outside the crit- 
ical line of fluxoid state [-l<^>o] but well within that of 
[O^o]- If the system is in state [-l^^^o], it will switch, pro- 
ducing a voltage which is detected. If it is in state [O^o], 
it will not switch. We achieve a one-shot discrimination 
between flux states of 100% with this readout. The pro- 
cess is completed by bringing to zero and waiting 100 
/is for the heat generated by a switching event to dissi- 
pate before repeating. Multiple repetitions, at a rate of 
about 5 kHz, yield Pesc- 

The overlapping Pesc = 50% critical lines seperating 
the [— l^o] and [O^o] fluxoid states are plotted in Fig.|2]3. 
Each ends in a cusp at the extreme value of flux where 



the corresponding fluxoid state is stable. These cusps 
occur at a non-zero current bias = ±alc = ±81 nA 
due to the critical current asymmetry a. The horizontal 
separation of the cusps scales precisely with 1/b. Our 
generalized MQT theory is accurately able to reproduce 
the measured data of flg|2]3. Of the parameters that go 
into this theory, b and a are treated as free parameters 
in this flt, Ic and rj are determined by the flt in Fig. [2^ 
and C is determined by a flt to spectroscopic data. 

Along the critical line of a given fluxoid state, for 
above or below the value /^^^p^ the potential is tilted to 
the right or to the left, and escape occurs preferentially in 
that direction. At /^^^p^ the camel potential is symmet- 
ric around the minima (a = 0), the two potential barrier 
heights are equal, and escape occurs with equal probabil- 
ity in either direction. The cusps in Fig. ^ correspond 
therefore to a double-path escape. 

The width of the escape process contains additional 
information about the dependence of the potential on 
the bias parameters, and on fluctuations in the bias pa- 
rameters [M]. In Fig. |2]3, we plot the width A$ = 
1^80% ~ ^20% 1 7 as a function of J^. This plot peaks around 
except that at this point there is a sharp dip (see 
insert). This behavior is explained by double-path MQT 
escape if we include low frequency current fluctuations. 
In this circuit thermal fluctuations are expected in J^, 
which we estimate to be on the order of 10 nA RMS by 
the equipartition theorem ^kT = ^P/rmS' where k is 
Boltzmann's constant, T :^ 40 mK is the circuit temper- 
ature, and P :^ 10 nH is the series isolating inductance. 
Because of this noise, the escape probability is averaged: 
(Pesc(^6, ^ext))- The angle brackets represent a convo- 
lution with the probability distribution of which we 
assume to be Gaussian with standard deviation /rms- 
As shown in Fig. |2]3, the addition of /rms = 9 nA is 
accurately able to explain both the increase in the over- 
all width, and the presence of a distinctive dip at /^^^^ 
which is a result of symmetry in escape direction. The 
presence of the dip and our ability to reproduce it with 
MQT theory is a striking conflrmation of double path 
escape and low frequency lb fluctuations in our sample. 

In Fig. |3^ and b we investigate the operation of a 
camel-back potential phase qubit corresponding to the 
two lowest levels |0) and |1) of the anharmonic central 
well related to the [O^o] flux state (see flg [l]^). For 
these measurements we use the same procedure as for the 
ground-state escape measurements except that before the 
nanosecond measurement pulse, an adiabatic flux pulse 
brings the system to the working point flux where 
a MW pulse is applied to the fast current line. At the 
working point, the barriers are high enough that Pesc is 
negligible. Immediately following the MW pulse we ap- 
ply a 5 ns pulse which projects the qubit state onto the 
flux state of the SQUID. This is possible because Pesc de- 
pends exponentially on the excitation level of the qubit. 
The amplitude of the measuring flux pulse is tuned such 
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FIG. 3: (a) Pesc versus h and MW frequency. Dark and 
bright grayscale correspond to high and small Pesc • (b) Width 
of the resonance on a semi-log scale. The dashed line is the 
predicted contribution due to 9 nA RMS low-frequency cur- 
rent noise. The dotted line is for 40 /i$o RMS low-frequency 
flux noise. The sum of these two contributions, the solid 
line, accurately reproduces the data (symbols). Rabi (c) and 
Ramsey (d) oscillations at the optimal line at h = —71 nA, 
^^^t/^o = -0.468 for the [-l$o] flux state. 



that escape will occur with high probability if the qubit 
is excited, and low probability if it is not. The measure- 
ment pulse transfers the quantum states |0) and |1) of the 
qubit to the classical fiuxoid states [O^o] and [— l^o] of 
the SQUID. Readout of the fiuxoid state, which is itself 
stable, reveals the projected qubit state, and repetition 
yields Pesc- 

Pesc was measured as a function of MW frequency u 
and lb. Because the MW pulse duration, 800 ns, is much 
longer that the relaxation time Ti :^ 100 ns, the system 
reaches a steady state. A peak in Pesc appears when 
1/ matches the qubit transition frequency vqi. Fig. [3^ 
shows Vol as function of I^. It reaches a maximum at 
/^^(<l>ext) which corresponds to the camel potential sym- 
metric point. Note that this optimal point is a function 
of flux, and is terminated by the cusp at the critical line. 
This data was taken at $ext = -0.503<l>o for the [0$o] 
fiuxoid state. Apparent in this spectroscopic image are 
avoided level crossings with what are likely microscopic 
two-level fluctuators, as first observed by Ref.[3 . We ob- 
serve on average 20 crossings per GHz. In Fig. Isb, the 



spectroscopic width of the uqi transition Ai^oi is plotted 
as a function of J^. A sharp minimum is observed at 
lb = 108 nA, corresponding to the flat maximum in uqi. 

We find that we can accurately model Ai^oi(A) with 
a combination of low-frequency current and flux fluctua- 
tions. Because vqi depends on the bias parameters, fluc- 
tuations cause uoi to vary from repetition to repetition, 
smearing out the observed resonance. Assuming a Gaus- 
sian fluctuation distribution, the predicted variance in 
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rent fluctuations alone, and (^f^)^ = 
for flux fluctuations alone. Here Ai^/ has been expanded 
to second order in /rms since ^ is zero at the opti- 
mal line. In Fig. [sJd, the predicted Az^/ is plotted as a 
dashed line for /rms = 9 nA, precisely the same current 
fluctuation amplitude used in Fig. The dotted line 
plots for <^>RMS = 40 jJi^o. The solid line is the com- 
bined prediction Az^ = ^J~Kl7f^^~Kl^. The dashed line 
is obscured behind the solid line except in a small region 
around the optimal current. This plot vividly demon- 
strates the idea of the optimal line: the effects of cur- 
rent bias fluctuations, which accurately account for the 
spectral width away from the optimal line, are rendered 
negligible on the optimal line. The residual spectroscopic 
width, about 10 MHz, can be explained by a flux noise 
of 40 Q RMS. Since the decoherence time T2 scales in- 
versely with Az^oi7 this optimal line is also optimal for 
qubit operations. Along this line Rabi and Ramsey oscil- 
lations (Fig. ^ and d) were measured giving coherence 
times of TRabi = 67 ns and pRamsey = 18 ns for this cur- 
rent sample. The anharmonicity is large enough and the 
applied power small enough that excitation beyond the 
first excited state is negligible, as we have verified by the 
linearity of Rabi frequency versus power. The system 
is confined to its lowest two levels and can therefore be 
considered a qubit. 

In conclusion, we have studied the quantum dynamics 
of a novel quadratic-quartic "camel" potential created in 
a dc SQUID circuit with 4 0, <l>ext — 0.5<l>o. Ground 
state escape exhibits critical line cusps and a dip in the 
escape width versus bias-current. We explain these two 
effects with a generalized double-path MQT escape the- 
ory. Moreover due to the particular potential symme- 
try, the quantum dynamics is insensitive in first order 
to current fluctuations along an optimal line /^^(^ext)- 
Along this line, the dc SQUID can be used as a phase 
qubit whose main decoherence sources are residual flux 
noise and microscopic two- level fluctuators. Future op- 
timization and exploitation of the unique properties of 
this system will aid in the understanding of decoherence 
mechanisms in quantum circuits and has the potential to 
yield a competitive phase qubit. 
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